Bloch Oscillation under a Bichromatic Laser: Quasi-Miniband Formation, Collapse, 
and Dynamical Derealization and Localization 
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A novel DC and AC driving configuration is proposed for 
semiconductor superlattices, in which the THz AC driving 
is provided by an intense bichromatic cw laser. The two 
components of the laser, usually in the visible light range, 
are near but not exactly resonant with interband Wannier- 
Stark transitions, and their frequency difference equals the 
Wannier-Stark ladder spacing. Multi-photon processes with 
the intermediate states in the conduction (valence) band cause 
dynamical derealization and localization of valence (conduc- 
tion) electrons, and the corresponding formation and collapse 
of the quasi-minibands. 

PACS numbers: 73.20.Dx, 42. 50. Hz, 72.20.Ht 



Semiconductor superlattices (SSL's) JD provide a stage 
for low-dimensional physics and for a variety of ultrafast 
optoelectronic devices. A typical example for such sys- 
tems is the well-known Bloch oscillation (BO) M and its 
frequency counterpart, the Wannier-Stark ladder (WSL) 
||. Mainly owing to their long artificial periods, SSL's 
have been the first systems in which the WSL and BO 
were observed B, verifying the early predictions [00] 
based on the band theory M. Since the BO is accom- 
panied by a time-dependent dipole momentum, biased 
SSL's are expected to generate tunable coherent electro- 
magnetic radiation in the THz range H , which may find 
extensive applications, e.g., in medical diagnoses. How- 
ever, the THz emission from Bloch oscillators has turned 
out to be quite weak, and a strong pulsed laser is re- 
quired for exciting the wave packets J6|, which limits its 
potential applications. 

A THz field has been used to drive the BO, as a re- 
sult the inverse Bloch oscillators can emit harmonics of 
the driving field JTJJ^]- As already reported, the biased 
SSL's subject to AC-fields exhibit a wealth of interesting 
physical effects p~pd| . The Hamiltonian of such systems 
is time-periodic, thus the eigenstates are the Floquet 
states associated with Floquet indices, i.e. quasiener- 
gies, forming the minibands p-pd[. With variation of 
the THz field strength, the quasi-miniband width oscil- 
lates, and will collapse at certain field strengths [p|-pd|. 
Correspondingly, the Bloch electrons are dynamically de- 
localized and localized |p|, p~0|p~^ [T^] . Experiments about 
these effects, however, have been quite scarce Rpf, be- 
cause it is not easy to obtain an intense THz source such 
as the free electron laser, which has moreover to be effec- 
tively coupled to the SSL's. This also prevents utilization 



of the THz driven SSL's in practical devices. 

In this Letter, we propose an alternative driving config- 
uration. As schematically illustrated in Fig. [y, a biased 
SSL is driven by an intense bichromatic cw laser with its 
two components at the frequencies u>i and 102 , and whose 
average $7 = [u\ + u>2)/2 is nearly resonant with E g , the 
frequency of the ground-state excitonic Wannier-Stark 
transition in the same quantum well, and thus is usually 
in the visible light range. The frequency difference be- 
tween the two components 2uj = \w\ — W2I is adjusted 
to ujboi th° BO frequency or WSL spacing, which is in 
the THz range. Suppose the two components have equal 
weight, then the optical field can be expressed as 



E(£) = 2E cos(w£)cos(rj£). 



(1) 



Here we have assumed the sample to be optically thin and 
neglected the spatial dependence of the laser. As will be 
shown later, this SSL can be equivalently described by 
a time-periodic Hamiltonian with the period 2tt/lu, and 
in which the dynamical effects are similar to those in the 
usual DC-biased and AC driven SSL's g-Q. 

At this point, it should be pointed out that although, 
atoms excited by two optical fields have been extensively 
studied both theoretically and experimentally pa], and 
the AC Stark effect in semiconductors driven by an in- 
tense monochromatic laser is also well understood [pi, 
to the best of our knowledge, there has been no report 
on Bloch oscillations driven by a bichromatic laser. 

In order to keep the following treatment as simple as 
possible, we have made several simplifications. First, we 
consider only the ground-state exciton for each Wannier- 
Stark transition. This is justified because the laser fre- 
quencies can be chosen well below the continuum transi- 
tions in the same quantum well and also by the fact that 
the ground-state exciton possesses by far the strongest 
component of oscillator strength. Secondly, we neglect 
the unequal WSL spacing effect |||l7| and Fano inter- 
ference between the continua and the embedded discrete 
states jl8| induced by the Coulomb interaction. Usually, 
the excitonic coupling affects significantly the dynamical 
derealization or localization |13Jl4], this however will 
be taken into account in our later work. Thirdly, the 
Zener tunneling and band mixing effects are not included. 
And lastly, the laser is assumed detuned for transitions 
other than those between the highest valence and the low- 
est conduction minibands. All these assumptions will in 
some way affect our results. Even so, the principle as well 
as main conclusions will not be qualitatively modified. 



With these simplifications, the Hamiltonian can be 
mapped into a one-dimensional model with two-band as 
well as nearest-neighbor coupling. With a well-justified 
rotating wave approximation for the optical coupling and 
the unitary transformation 

S = exp[irtt^(al n a cn - al n a vn )/2], 

I! 

the Hamiltanian in the rotating frame is 
H = 22 (™ BO + y ) a c« a cn - -^-a\ n {a cn+ i + a c „_i) 

n 

+ (nuj B o - y) a La v „ + — -al„(a vn+ i + a v „_i) 



-X cos(ujt) (al n a vn + H.c. ) , 



(2) 



where a\ n (aj. n ) stands for the operator creating a con- 
duction (valence) electron at the nth site, A c (A v ) is the 
conduction (valence) miniband width, eo = E g — Q is the 
detuning, and \ = cIEq is the interband optical transition 
strength (d is the interband dipole matrix element along 
the E direction). Compared to the usual THz driving 
WSL system, the driving force in the present configura- 
tion is essentially an interband rather than an intraband 
process. 

The time-periodic system described by Eq. (g) has no 
stationary solution, and the eigenstates are time-periodic 
Floquet states \q, m) = exp(imujt)\q) with period T = 
2tt/lu, which satisfy 



(H -id t )\q,m) 



-q,m 



q,m), 



where e q , m = e q + muj is the Floquet index, i.e. 
quasienergy. The Floquet states and indices can be ob- 
tained by diagonalizing the secular equation 

U(T,0)\q)=exp{-ie q T)\q), 

where the propagator 



U{T,0) = fexp[-i / H(t)dt] 



The Hamiltonian (|2|) can be decomposed into a sum 
of mutually commutative 2x2 matrices in the accel- 
erating quasimomentum representation {|a, k)\a —c, v} 
with the transformation \a, k) — \a, k — LUsot) (hereafter 
the superlattice period is assumed to be unity). That is 
H = Ek H k, with 

H * = 1] £ a k a ik a *k ~ XcosM)(o^Ovfc + H.c), (3) 



where £ c j. 



(eo — A c cosfc)/2 and e v ^, = (— eq + 
A v cos k)/2 are the conduction and valence miniband dis- 
persions, respectively. Thus each quasimomentum k cor- 
responds to two sets of Floquet states \±,m)k with the 
quasienergies £k±, m - 



Before we present our numerical results, the quasi- 
minibands are analytically evaluated in the weak and 
strong optical coupling limits, so as to give some physical 



insight. From now on, we will set E„ 



Q. 



When x = 0, i.e. there is no optical field, Hk is al- 
ready diagonalized, and the quasienergies produce the 
well-known WSL §. 

When x is small, the optical coupling 

Him = -xcos(wi) (a^a vS + H.c.) 
can be treated as a perturbation to 

#o,fc = e c fe a lfe«cfe + £ v ~k a l~ k a v~k- 

Similar to the steady case, a perturbation theory can be 
developed, in which the matrix element for time-periodic 
functions is defined as 

((q,m\H\q ,m)) = -I (q,m\H\q ,m)dt 

in the Hubert subspace |H]]. In terms of unperturbed 
Floquet states 

|±) fc = exp[±iA c/v /(2w BO ) sin(cj BO t - fc)]|c/v, fc) 

with the quasienergies e k ± = 0, the optical coupling is 
expressed as 

k {{+,m\H^ k \+,m')) k = k ({-,m\H ltk \-,m'))k = 0, (4) 
k {{-, m\Hi, k \+,m'))k =k ((+,m'\H 1<k \-,m))* k 



9 Z-^z n,- 



M 



A 



n.± 



2ujbo 



)e- 



(5) 



where A = A c + A v is the combined miniband width, n 
is an integer number, and J n (x) is the nth order Bessel 
function of the first kind. According to Eq. (Q) and (||), 
the quasienergy spectrum depends on only A rather than 
the ratio A c /A v . 

From the symmetry of the matrix elements, the 
quasienergies e k ± satisfy the relation e k ± = — £fc T = 
£-K-k^ = £ fc-7TT = £ -k±, so each Floquet state is 4-fold 
degenerate, and only the quasienergies for k € [0, n/2] 
need to be calculated to obtain the full spectrum infor- 
mation. 

Since the first-order modification for the quasienergies 
vanishes (Eq. (H)), the second-order perturbation gives 



£fc± = TX 2 cosk sr ^ 



+ oo j i A \ j f_A 
Jn y2ujBo' n ~ 1( - 

(2n - 1)lo 



2ujbo ' 



n=\ 



(G) 



which forms a miniband. It can be recognized at once 
that, except for a prefactor, the quasi-miniband has 
a dispersion relation identical to that of the original 
minibands. The discrete WSL is broadened, indicat- 
ing dynamical derealization of the originally localized 
Wannier-Stark states. 



An intuitive picture favors our understanding of the 
dynamical derealization. As depicted in Fig. [j] (b), 
a valence band electron in a Wannier-Stark state can 
hop to its neighbor states by first absorbing a photon 
from one component of the bichromatic laser and then 
emitting a photon to the other component, though both 
components are off-resonant with the interband transi- 
tions. The square law dependence of the quasi-miniband 
width on the optical coupling strength \ results from this 
second-order process, in which a conduction band elec- 
tron plays the role of an intermediate state. Clearly, the 

Bessel functions Jni^^) and ^-1(2^) in Ec l- ® 
reflect the two virtual transitions between the intermedi- 
ate state and the initial and final Wannier-Stark states. 
A conduction electron can similarly be dclocalized. 

With increasing laser intensity, multi-photon processes 
become more and more important, and the perturbation 
treatment for the optical coupling will fail eventually. 
Thus, when x ^ A, the Hq^ should be regarded as the 
perturbation to H\.k |10|| . Since eigenstates of H\^ are 

\±)k = exp[±ix/wsin(ui)]%/2/2(|c,fc) ± |v,fc», 

the first-order degenerate perturbation gives the quasi- 
miniband dispersion as 



£fc± = ± — J2{ — )cosfc. 



(7) 



In this strong coupling limit, the quasi-miniband width 
oscillates according to a second order Bessel function 
with the argument of 2\ over lo, and the quasi-miniband 
collapses at such laser intensities that 2x/w is a root of 
the Bessel function, which is very similar to the case of 
intra-miniband driving by a THz field. The Taylor ex- 
pansion of the Bessel function in Eq. (f7|) contains only 
even powers of x, indicating that each virtual photon ab- 
sorption is necessarily followed by a virtual photon emis- 
sion, and no process with an odd number of photons is 
involved. This again confirms our intuitive picture in- 
voked to understand the dynamical derealization. 

The formation and collapse of the quasi-miniband can 
be clearly seen from Fig. @, in which the quasiener- 
gies are numerically calculated as functions of the op- 
tical coupling strength for 21 k points taken evenly from 
to 7r/2. For comparison, the perturbation results for 
k = are also plotted. In the weak optical coupling 
limit (2%/w < 2), the second order perturbation of Hi t k 
(dashed lines) gives results well consistent with the nu- 
merical ones. In the strong coupling limit, the first order 
perturbation of i?o,fc is also a quite good approximation, 
especially for smaller A, where the requirement for the 
perturbation treatment x 3> A is better met. 

Comparing the numerical results with the perturba- 
tion calculations, two points need mentioning. First, the 
quasi-miniband collapses at 2x/w «4.16, 7.84, etc. for 
A = 33 meV, and at 2 X /lo «4.96, 8.32, etc. for A = 16.5 



meV, which correspond to laser intensities a little weaker 
than the roots of the Bessel function. We think it comes 
from the underestimated localization effect of the DC 
field by treating Hq^ as a perturbation. Secondly, the 
numerically calculated band collapse is not as complete 
as that in perturbation calculations, because of the inter- 
action between different sidebands of the Floquet states 
of Hi,fc. 

To show the dynamical derealization and localization 
of a bichromatic laser driven electron, we have calculated 
the evolution of an initially localized electron by integrat- 
ing the Shrodinger equation. The wave packet width W 
defined to characterize the localization of an electron is 



W 2 



Z^ 



(n - r ) 2 {(al n a cn ) + (4„a vn )), 



where 



r = 2^ n 



<(4 



^))■ 



When t = 0, the electron is put at the valence site n = 
]19[ . With such an initial condition, the wave packet 
center tq is found to be always zero for any laser intensity. 

In Fig. ||, W is shown as a function of time within the 
first ten periods for three laser intensities. Without the 
optical field (Fig. g (a)), the wave packet width W oscil- 
lates sinusoidally with time with a static center-of-mass, 
which is just the breathing mode of the BO pQ| and the 
well-known Wannier-Stark localization. The two non- 
zero laser intensities correspond to approximately max- 
imum quasi-miniband width and band collapse, respec- 
tively. When the quasi-miniband width is a maximum 
(2x/u> — 3.04, Fig. (b)), W increases linearly with 
time with a slight oscillation superimposed, which clearly 
evidences the dynamical derealization of the initially 
localized electron. When the quasi-miniband collapses 
(2x/c<j = 4.96, Fig. g (c)), W oscillates with time signifi- 
cantly and quasi-periodically, whose minima are slightly 
shifted upward due to imperfect band collapse. It is ob- 
vious that the dynamical localization is realized. 

Note that the dynamical localization length is larger 
than the Wannier-Stark localization length by a factor 
of about A/A v , which suggests a similar band- width de- 
pendence for both kinds of localization. This behavior 
is also exhibited for dynamical derealization, which, as 
shown in Fig. (b), compares very well to the diffusion 
of a free Bloch electron provided identical bandwidths 
(«3.9 meV) are adopted. 

The population of the conduction miniband f c = 
En( cn a cn) is also plotted for 2x/w = 3.04 (the dot- 
ted line). / c oscillates nearly periodically with time due 
to virtual excitation. When a valence band electron is 
excited into the conduction states, it can hop to neigh- 
bor sites, thus W ascends by a step, which confirms once 
more our intuitive picture for the dynamical derealiza- 
tion process. 



In summary, we have predicted the quasi-miniband for- 
mation and collapse, and correspondingly the dynami- 
cal derealization and localization for an electron in a 
SSL driven by DC and AC fields in a novel configura- 
tion, in which the THz AC driving is provided by an 
intense bichromatic visible laser with the frequency dif- 
ference of the two components equal to the WSL spacing. 
Our findings suggest that the experiments on the dynam- 
ical effects can be carried out without expensive THz 
sources like free electron lasers. Suppose lobo = 2o> = 10 
meV, an optical coupling of 12.4 meV is enough to re- 
alize the dynamical localization, which should be acces- 
sible already in many laboratories. It should be noted 
that the two components of the bichromatic laser can be 
replaced by two independent lasers with small phase fluc- 
tuation. The possible THz emission from so driven Bloch 
oscillators, requiring neither an intense THz source nor 
pulse optical excitation, may find broad applications. In 
addition, this interband driving configuration introduces 
several interesting physical problems such as the exci- 
tonic effect and excitation-induced many body correla- 
tion, which are absent in the usual intraminiband THz 
driving configuration. 
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FIG. 3. Wave packet width (solid lines) as functions of time 
within the first ten periods for 2x/w = 0, 3.04, and 4.96 in 
(a), (b), and (c), respectively. Other parameters are the same 
as in Fig. (b). In (b), the dashed line stands for wave packet 
width of a free Bloch electron which is initially localized, and 
the dotted line is the conduction electron population. 
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